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The surface Fermi arc states in Z 2 Weyl semimetals AsBi (A = Na, K, Rb) are studied by 
employing a continuum low-energy effective model. It is shown that the surface Fermi arc states 
can be classified with respect to the ud-parity symmetry. Because of the symmetry, the arcs come 
in mirror symmetric pairs. The effects of symmetry breaking terms on the structure of the Fermi 
arcs are also studied. Among other results, we find at least two qualitatively different types of the 
surface Fermi arcs. The arcs of the first type link disconnected sheets of the bulk Fermi surface, 
while arcs of the second type link different points of the same bulk Fermi surface sheet. 

PACS numbers: 73.20.At, 71.10.-w, 03.65.Vf 


I. INTRODUCTION 


3D Dirac semimetals are 3D analogs of graphene [ij. Their conduction and valence bands touch only at discrete 
(Dirac) points in the Brillouin zone with the electron states described by the 3D massless Dirac equation. Each Dirac 
point in momentum space is composed of two superimposed Weyl nodes of opposite chirality. Such points are usually 
obtained by fine tuning of certain physical parameters (e.g., the spin-orbit coupling strength or chemical composition) 
and are difhcult to control. Additionally, they are often unstable with respect to mixing of Weyl modes and opening 
a gap. 

An important idea was proposed in Refs. (2|, [Sj, where it was shown that an appropriate crystal symmetry can 
protect and stabilize the gapless 3D Dirac points. Indeed, if a pair of crossing bands belong to different irreducible 
representations of the discrete (rotational) crystal symmetry and if this symmetry is not broken dynamically, then the 
mass term for the corresponding Dirac fermions will be prohibited. The ab initio calculations in Ref. Q showed that 
/3-cristobalite Bi02 exhibits three Dirac points at the Eermi level. Unfortunately, this material is metastable. By using 
the first-principles calculations and an effective model analysis, the compounds AsBi (A=Na, K, Rb) and Cd 3 As 2 
were identified in Refs. ii as possible 3D Dirac semimetals protected by crystal symmetry. Giant diamagnetism, 
linear quantum magnetoresistance, and the quantum spin Hall effect are expected in these materials. Furthermore, 
various topologically distinct phases can be realized in these compounds by breaking the time-reversal and inversion 
symmetries. By using angle-resolved photoemission spectroscopy, the Dirac semimetal band structure was indeed 
observed in Gd 3 As 2 and Na 3 Bi opening the path toward experimental investigations of the properties of 3D 
Dirac semimetals. 

Weyl semimetals is another group of materials that is closely related to 3D Dirac semimetals and have already 
attracted a lot of theoretical interest (for reviews, see Refs. ). They are characterized by topologically non¬ 

trivial Weyl nodes in reciprocal space. Weyl nodes are monopoles of the Berry flux and, therefore, can appear 
or annihilate only in pairs. Weyl semimetals were proposed to be realized in pyrochlore iridates [l^ . topolo gica l 
heterostructures [l^ , magnetically doped topological insulators 0 , and nonmagnetic materials such as TaAs [l5lll6j| . 
Recently, first experimental studies of Weyl semimetal candidate TaAs were reported in Refs. [TtI - I^ . The authors 
observed unusual transport properties and surface states that are characteristic of the Weyl semimetal phase. Another 
interesting realization of the Weyl points in the context of photonic crystals has been recently reported in Ref. 211. 

Since the magnetic field breaks the time-reversal symmetry, a Dirac (semi-)metal in a magnetic field may transform 
into a Weyl one with Weyl nodes separated in momentum space by a nonzero chiral shift [ 2 ^ . Experimentally, the 
transition from a Dirac metal to a Weyl one in a magnetic held might have been observed in Bii_xSbx for x pz 0.03 
23|. In moderately strong magnetic helds, a negative magnetoresistivity is observed and interpreted as a hngerprint 
24 -[ 2 ^ of a Weyl/Dirac metal phase. 

The surface Fermi arcs [1^ [27H^ , which connect Weyl nodes of opposite chirality, are related to the non-trivial 
topology of Weyl semimetals. In equilibrium, the presence of such surface states ensures that the chemical potentials 
at different Weyl points are identical Although Fermi arcs always connect Weyl nodes of opposite chirality, 
their shapes depend on the boundary conditions and, as shown in Ref. (30j |. Fermi arcs of an arbitrary form can 
be engineered. The Fermi arcs on the opposite surfaces of a semimetal sample together with the Fermi surfaces of 
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bulk states form a closed Fermi surface. In an external magnetic field, the nontrivial structure of the corresponding 
Fermi surface gives rise to closed magnetic orbits involving the surface Fermi arcs (3ll | . These orbits produce periodic 
quantum oscillations of the density of states in a magnetic field leading to unconventional Fermiology of surface states. 
It was argued in Ref. [IH that the interaction effects can change the separation between Weyl nodes in momentum 
space and the length of the Fermi arcs in the reciprocal space and, thus, affect these magnetic orbits. As a result, 
we found that the period of oscillations of the density of states related to closed magnetic orbits involving Fermi arcs 
has a non-trivial dependence on the orientation of the magnetic field projection in the plane of the semimetal surface 
[s^. If experimentally observed, such a dependence would provide an important clue to the effects of interactions in 
Weyl semimetals. 

Normally, one would not expect any surface Fermi arcs in 3D Dirac semimetals because the Dirac point has 
no topological charge and the associated Berry flux vanishes. In Refs. Si , however, it was shown that the 3D 
Dirac semimetals AsBi (A=Na, K, Rb) and Cd 3 As 2 possess non-trivial surface Fermi arcs. This finding suggests a 
topologically nontrivial nature of the corresponding Dirac materials. Recently we showed that this is indeed the 
case for Dirac semimetals AaBi (A = Na, K,Rb). The physical reason for their nontrivial topological properties is 
connected with a discrete symmetry of the low-energy effective Hamiltonian. The symmetry classification allows one 
to split all electron states into two separate sectors, each describing a Weyl semimetal with a pair of Weyl nodes 
and broken time-reversal symmetry. The time-reversal symmetry is preserved in the complete theory because its 
transformation interchanges states from the two different sectors. The nontrivial topological structure of each sector 
was supported by explicit calculations of the Berry curvature, which revealed a pair of monopoles of the Berry flux at 
the positions of Weyl nodes in each of the two sectors of these semimetals . In essence, these results demonstrated 
that Dirac semimetals AaBi (A = Na, K,Rb) are, in fact, Z 2 Weyl semimetals. 

In Refs. ai, the surface Fermi arcs in 3D Dirac semimetals were obtained in a tight-binding model by using an 
iterative method that produces the surface Green’s function of the semi-infinite system [33| . The imaginary part of 
the surface Green’s function makes possible to determine the local density of states at the surface. While such a 
technique is very powerful, it is essentially a “black box”. In contrast, in the present paper, we study analytically the 
surface Fermi arc states by employing the continuum low-energy effective model with appropriate boundary conditions 
at the surface. We hope that such a consideration will provide a deeper understanding of the physical properties and 
characteristics of the surface Fermi arcs, as well as shed more light on the nontrivial topological properties of the 
AaBi compounds. 

The paper is organized as follows. In Sec. m we introduce the low-energy effective model and discuss its symmetries. 
The recently revealed Z 2 Weyl semimetal structure of AsBi (A = Na, K, Rb) is emphasized. In order to clarify the 
origin and the structure of the surface Fermi arcs, we study in Sec. Illll the corresponding states in a simplified model 
that contains a single Weyl semimetal sector. In Sec. lIVi we present the rigorous analysis of the surface Fermi 
arc states in a realistic low-energy model of semimetals AsBi (A = Na, K, Rb). The effects of several possible 
symmetry breaking terms on the structure of the surface Fermi arc states are investigated in Sec. [V] The discussion 
and the summary of the main results are given in Sec. IVII Technical details regarding the symmetry properties and 
classification of the Fermi arc states are presented in Appendices O and |B] 

For convenience, throughout the paper, we set h = 1 and c = 1. 


II. MODEL 

A. Low-energy effective Hamiltonian 

The low-energy Hamiltonian derived in Ref. Q for AaBi (A = Na, K, Rb) has the form 

H(k) = eo(k) + H4x4, 

where eo(k) = Cq + Cik^ + C2{k'^ + ky) and 

/ M(k) Ak+ 0 B*{k) \ 

_ Ak_ -M(k) H*(k) 0 1 

If4x4 - Q 

V B{k) 0 -Ak+ -M(k) J 


( 1 ) 


( 2 ) 


While the diagonal elements of ^ 4 x 4 are given in terms of a single function, M(k) = Mq — Mik^ — M2{k‘^ -I- ky), the 
off-diagonal elements are determined by functions Ak± and i?(k) = ak^k'^, where k± = kx A: iky. 
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By fitting the energy spectrum of the effective Hamiltonian with the ab initio calculations, the numerical values of 
parameters in the effective model were determined in Ref. They are 


Co = -0.06382 eV, 

Cl = 8.7536 eVA^ 

C 2 = -8.4008 eVA^ 


Mo = -0.08686 eV, 

Ml = -10.6424 eVA^ 

M 2 = -10.3610 eVA^ 

(3) 

A = 2.4598 eVA, 

a = 5.448 A, 

c = 9.655 A, 



where we also included the lattice constants a and c. Since no specific value for a was quoted in Ref. , we will treat 
it as a free parameter below. 

The energy eigenvalues of the low-energy Hamiltonian o are given by the following explicit expression: 

£;(k) = eo(k) ± VM2(k)+H2fc+fc_ + |H(k)|2. (4) 

It is easy to check that the term with the square root vanishes at the two Dirac points, k^ = (0,0,±i/m), where 
■y/m = With the choice of the low-energy parameters in Eq. ([3]), we find that y/m ^ 0.09034 A The 

function B(k) plays the role of a momentum dependent mass (gap) function that vanishes at the Dirac points. 

It is instructive to show that linearizing M(k) in the vicinity of the Dirac points k^, Hamiltonian ([2]) takes the 
form of a 3D massive Dirac Hamiltonian. In the vicinity of k))", expanding M(k) to the linear order in 5k = k — k))", 
we obtain 


rjlin _ 

-^4x4 — 


f - kyffy - k^az) i?*(k)cra, 


i?(k)c 


—A k • cr 


(5) 


where cr are Pauli matrices and k = (kx, ky, 25kz\/M qMi/A). Furthermore, by performing the unitary transformation, 
= Ux where Ux = diag{ax, I 2 ) and I 2 is the 2x2 unit matrix, we find that the Hamiltonian takes the 

standard form of the Dirac Hamiltonian in the chiral representation, 


ulin _ 

-^4x4 ~ 


Hk-cr B*(k) \ 

H(k) —A k ■ cr J 


( 6 ) 


Taking into account that the mass term B(k) vanishes at the Dirac point, we conclude that the upper and lower 
2x2 blocks describe quasiparticle states of opposite chiralities. Also, since the leading order nonzero corrections to 
the mass function are quadratic in momentum, the chirality remains a good quantum number in a sufficiently small 
vicinity of the Dirac point. Hamiltonian (jS]), describing two subsets of the opposite chirality states near a single Dirac 
point, does not appear to have any interesting topological properties. Also, by itself, it is unlikely to give rise to any 
Fermi arcs states. It is easy to check, however, that Hamiltonian Q linearized near k))" has a similar structure and 
describes two additional subsets of the opposite chirality states. As we argue below, the superposition of the two 
sectors of the theory is nontrivial and gives rise to an interesting topological structure [s^. 


B. Symmetries 


Let us briefly review the symmetry properties of the low-energy Hamiltonian following Ref. [^. We start by 
pointing out that, as expected, the Hamiltonian © is invariant under the time-reversal and inversion symmetries, 


i.e., 


0iJ_k0-^ = Bk, 
FB^kP-^ = iJk, 


(time-reversal symmetry) 
(inversion symmetry) 


(7) 

( 8 ) 


where Q = TK (K is complex conjugation) and 


/ 

0 

0 

1 


/I 

0 

0 



0 

0 

0 

1 1 p 

0 

-1 

0 

0 1 


-1 

0 

0 

0 ’ 

0 

0 

1 

0 

V 

0 

-1 

0 

0/ 


0 

0 

-1/ 


(9) 


Having both, the time-reversal and inversion symmetries, suggests that the corresponding compounds are not Weyl 
semimetals. This is not the whole story, however. 
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As shown in Ref. [s^, the low-energy Hamiltonian in Eq. o possesses a new discrete symmetry, the so-called 
up-down parity (ud-parity), that protects its topological nature. In order to understand the corresponding symmetry, 
it is instructive to start from the approximate Hamiltonian without the mass function R(k) (or, alternatively, a = 0). 
In this case, the 4x4 Hamiltonian takes a block diagonal form: 174x4(a = 0) = i7;^x2 ® ^ 2 x 2 - explicit form of 
the upper block is given by 

Mo — Ml 7^ — M2{kl + ky) A{kx + iky) 

A{k,-iky) -[M^-M^kl-M^ikl + kD] 

This block Hamiltonian defines a Weyl semimetal with two Weyl nodes located at kj. (The lower block H 2 y ^2 ^ 

similar form, except that k^ is replaced hy —k^-) It is well known [2^|3l| that such a Weyl semimetal has the surface 
Fermi arc connecting the Weyl nodes of opposite chirality at k^J" and k^j". Because of the sign difference, k^ —>■ —kx^ 
the chiralities of the states near the Weyl nodes at k(}^ are opposite for the upper and lower block Hamiltonians. Thus, 
the complete 4x4 block diagonal Hamiltonian 774 x 4 ( 0 ; = 0) describes two superimposed copies of Weyl semimetal 
with two pairs of overlapping nodes. Since the opposite chirality Weyl nodes coincide exactly in momentum space, 
they effectively give rise to a pair of Dirac points at k^. At the same time, because the opposite chirality nodes 
come from two different Weyl copies, they cannot annihilate and cannot form topologically trivial Dirac points. In 
fact, the corresponding approximate model describes a Z 2 Weyl semimetal [s^. The nontrivial topological properties, 
associated with the underlying Z 2 Weyl semimetal structure, ensure that the resulting Dirac semimetal possesses 
surface Fermi arcs. 

It is easy to show that the existence of the Z 2 Weyl semimetal structure in the absence of 77 (k) is connected with 
the continuous symmetry U-)-(l) x U_(l) of the approximate Hamiltonian 774 x 4 ( 0 ; = 0). This symmetry describes 
independent phase transformations of the spinors that correspond to the up- and down-block Hamiltonians, 77 ^^2 
77^^ 2 ) respectively. 

For 7?(k) ^ 0, the continuous symmetry U+(I) x U_(I) is broken down to its diagonal subgroup Uem(l) that 
describes the usual charge conservation. However, the low-energy Hamiltonian © with the momentum dependent 
mass function 77(k) = akzk\ possesses a ud-parity, dehned by the following transformation [ 3 ^ : 

t/77_fc^t/"^ = 77fc^, (ud-parity), (II) 




where matrix U has the following block diagonal form: U = diag(72, —I 2 ) and I 2 is the 2x2 unit matrix. For the 
Hamiltonian to be symmetric under the ud-parity, it is crucial that the mass function 77 (k) changes its sign when 
kz —>■ —kz [while the functions eo(k) and M(k) in the diagonal elements do not change their signs]. In the special 
case of a momentum independent mass function, such a discrete symmetry does not exist. 

As was argued in Ref. [33| , the existence of the noncommuting time-reversal and ud-parity symmetries implies that 
the AaBi semimetal is, in fact, a Z 2 Weyl semimetal. In such a semimetal, all quasiparticle states can be split into 
two separate groups, labeled by the eigenvalues x = of where Hfc^ is the operator that changes the 

sign of the z component of momentum, kz —>■ —kz- Effectively, each group of states defines a Weyl semimetal with a 
broken time-reversal symmetry. The corresponding symmetry is preserved in the complete theory, in which the two 
copies of Weyl semimetals are superimposed. 

The Z 2 Weyl semimetal structure of AaBi (A = Na, K, Rb) is also supported by the explicit calculation of the Berry 
connection and the Berry curvature in each Weyl sector described |33j . In particular, the corresponding results for 
the curvature in the momentum space reveal a clear dipole structure. It is natural, that each Weyl sector, described 
by quasiparticle states with a fixed eigenvalue of 7/^, should give rise to Fermi arcs connecting the pairs of Weyl nodes 
at k^ . Moreover, such arcs should be topologically protected and could not be removed by small perturbations of 
model parameters. 

In our discussion of Fermi arcs below, it will be also useful to take into account that there exists yet another discrete 
symmetry defined by the following transformation: 


where 


UH.U-^ = 77fc 


U = 


/O 0 1 0 
0 0 0 1 
10 0 0 
Vo 1 0 0 


( 12 ) 


(13) 


It is interesting to note that the product of the and Ullk^ transformations is also a symmetry 

of the low-energy Hamiltonian ©■ The symmetry THfe^H^^ is related to the time-reversal symmetry. This follows 
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from the fact that KIVky is also the symmetry of the low-energy Hamiltonian O- Together the operators 
and THfe^Hfe^ form a non-commutative discrete group. 

Hamiltonian o is rather complicated, therefore, the corresponding analytic calculations of its surface Fermi states 
are quite involved and not much revealing. Therefore, our general strategy in analyzing these states will be to start 
from a simplified model and then move forward to the realistic model by adding step-by-step the necessary missing 
pieces. 


III. SURFACE FERMI ARCS IN SIMPLIFIED 2x2 MODEL 


In order to get an insight into the structure of the surface Fermi arcs in the low-energy model described by 
Hamiltonian it is instructive to first study the surface Fermi arcs in a simplified 2x2 model, given by one of the 
diagonal blocks, e.g., Elq. ITOl) . (The solutions for the other block Hamiltonian, i?;7x2i can be obtained simply 

by changing kx —>■ —kx-) For completeness, we will also include the term eo(k) proportional to the unit matrix, which 
is present in the low-energy Hamiltonian. Thus, our model 2x2 Hamiltonian reads 


^2x2 — eo(k) + ^2x2 ~ Co(k) -I- 


Mo - _ M2{kl + kl) 

Aikx - iky) 


A (kx “t“ i ky ) 
[Mq — Mikl — M2{kl 


Before proceeding to the analysis, it is convenient to perform a unitary transformation, H 2 x 2 = 
Uy = -^ {I 2 + i'^y)- The transformed Hamiltonian has the following explicit form: 



U-^H2x2Uy, 


(14) 

where 


H 2 x 2 = eo(k) -k [7 [kl -m) - M 2 {kl + kl)] Ox - vkxd^ - vkyay, (15) 

where we introduced the notations similar to those in Ref. [ 2 ^: v = A and 7 = —Mi. 

To study the surface Fermi arcs, we will assume that the surface of a semimetal is at y = 0. The semimetal itself 
is in the upper y > 0 (lower y < 0) half-space when we describe the surface arc states on the bottom (top) surface. 
(Of course, in the absence of any effects that break the inversion symmetry ky —>■ —ky explicitly, the two cases will 
be related by a simple symmetry transformation.) Without loss of generality, we will concentrate primarily on the 
bottom surface states. The boundary condition on the semimetal surface will be imposed by replacing the parameter 
m with the — m on the vacuum side of the boundary and taking the limit m —> 00 [ 2 ^. From a physics viewpoint, 
such a replacement is the simplest way to prevent quasiparticle from escaping into the vacuum. 

Taking into account that the Fermi arc states should be localized at the y = 0 boundary, let us rewrite Hamiltonian 
dm) in the following form: 


fr _( Co + Cikl + C2{kl - dl) - vkx 'y {k^ - m) - M2(kl - dl) + vdy \ 
\1 {kl — m) — M2(fc^ — dl) — vdy Co + Cikl -I- C2(fc^ — 9 ^) -I- vkx j 

where, for the convenience of further derivations, we replaced ky = —idy. 


A. Simplified model with C 2 = M 2 = 0 


We will see in what follows that the presence of the terms with the second derivative with respect to y in Hamiltonian 
(fTSll leads to many technical complications and makes the analysis rather involved. Therefore, to set up the stage, in 
this subsection we start our analysis in an even more simplified model, described by Hamiltonian (HU) with C 2 and 
M 2 set to zero. Then, by introducing the two-component spinor 4' = , we see that the eigenvalue problem 

{H 2 x 2 — £')4' = 0 is equivalent to the following system of equations: 

{-vkx + Cikl + Co) tpi + [vdy + jkl - "fm{y)]'ip 2 = Eipi, (17) 

{vkx + Cikl + Co) iIj 2 + [-vdy + ykl --fm{y)] tpi = E'>p2- (18) 

Here m(y) = mOijj) — m0{—y), where 9{y) is the step function. Recall that, by assumption, the boundary condition 

at y = 0 is enforced by taking the limit m —>■ 00 on the vacuum side (y < 0). Formally, Eqs. (1171) and (fTSl) have the 

following surface state solutions: 


^ 1 ( 2 /) 


0 j’ 


^2(y) 




(19) 
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In the region occupied by the seiniinetal {y > 0), the solution 'l'i(y) is normalizable only for k'^ — m < 0, while the 
solution ^' 2 ( 2 /) is normalizable only for — m > 0. However, on the vacuum side (y < 0), only is normalizable. 

The dispersion relation for this normalizable surface state solution follows from Eq. (1171) . It is given by 

E = -vk^ + Cikl + Co- (20) 


By making use of this relation, we derive the equation for the bottom surface Fermi arc in the transverse k^k^ plane. 


k 


X — 


E-Cikl-Co 

V 


( 21 ) 


It is instructive to compare this surface Fermi arc with that in the model of Ref. [ 2 ^, where Ci = 0. While the surface 
Fermi arcs run between k^ = —\/m and kz = in both models, the arcs in the model of Ref. [ 2 ^ do not depend 
on the momentum k^- This is in contrast to the surface Fermi arc in Eq. (EU, for which kx is a quadratic function 
of kz- Thus, we see that the presence of the quadratic in kz term in the diagonal component of Hamiltonian (1161) 
produces a nonzero curvature of the surface Fermi arcs in momentum space. For illustration, several surface Fermi 
arcs for different values of the Fermi energy are shown in Fig. [TJ The arcs have parabolic shapes. The corresponding 
arcs in the model of Ref. [ 2 ^ would be given by straight lines. 



FIG. 1: (Color online) The bottom surface Fermi arcs for several different values of the Fermi energy in a simplified two- 
component model, described by Hamiltonian (ESJ with C 2 = M 2 = 0. The analytical form of the arcs is given in Eq. E3. 


Before concluding this section, let us note that the solution '^ 2 {y) in Eq. (fT^ describes Fermi arcs on the top 
surface. We find from Eq. (ITSl) that the corresponding dispersion relation is given hy E = vkx + Cikl + Cq. Let us 
also note in passing that there exists another set of the (top and bottom) Fermi arcs for the lower block Hamiltonian, 
^ 2 x 2 - The corresponding arcs are obtained from the solutions for the upper block Hamiltonian, iL^^ 2 i by making the 
replacement kx —>• —kx- 


B. The case with C 2 0 and M 2 yf 0 

Let us now consider the general case with C 2 0 and M 2 ^ 0. By noting that the Hamiltonian in Eq. dUl) contains 
second derivatives with respect to y, the eigenvalue problem ( 1 / 2 x 2 — E)'i> = 0 becomes more complicated. In the 
semimetal (j/ > 0 ), it is equivalent to the following system of coupled equations: 

[C2{kl - dy) - vkx + Cikl E Co] tpi + [-M2{kl - dy) + vdy E ykj - 7 m] ^’2 = Etpi, 

[C2{kl - dy) + vkx + Cikl + Co] tp2 + [--^ 2(^1 - d^) - vdy + ykl - ym\ = Eil)2. 


( 22 ) 

(23) 
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FIG. 2: (Color online) The bottom surface Fermi arcs (I25II for several values of the Fermi energy in a two-component model, 
described by Hamiltonian m with (72 0 and M 2 yf 0. 


On the vacuum side {y < 0), the corresponding set of equations has the same form, but with m replaced by —m. 
At the vacuum-semimetal interface (y = 0), the wave functions and their derivatives should satisfy the conditions of 
continuity, see Eqs. m through (IA4I) in Appendix lA II 

The key details of the derivation of the surface Fermi arc solutions are presented in Appendix [XTJ On the semimetal 
side, the spinor structure of the solution takes the following form: 

i=l k ® / 

where the explicit expressions for the exponents are given in Eq. (IA9I) . Note that the exponents take real values in 
the case of surface Fermi arc states. The condition of existence of nontrivial surface Fermi arc solutions is given by 

-C2{pI - kl) + Cikl + Cq - E - vkx _ -C2{pl - fcg) -b Cikl + Cp - E - vk^ 

-M 2 {pI - kl) - - m) + vpi -M 2 {p 2 - kl) - j{kl - m) + vp 2 

This equation defines the functional dependence kz{kx) for the possible surface Fermi arc states. A numerical study 
shows that nontrivial solutions exist only in a finite range of energies, i.e., —0.168 eV < A < 0.373 eV. Several 
solutions for different values of the energy are shown in Fig. [2] The results of the numerical analysis show that the 
following condition is satisfied: bi/ai = 62/02 = 0.5115 for all solutions. It is worth noting that the E = 0 surface 
Fermi arc in Fig. [5] appears to be almost identical to the corresponding arc, obtained by a very different method in 
Ref. 0, see Fig. 3c in that paper. 

So far, we considered the arc states only for one of the two-component block Hamiltonians, defined in Eq. (Hi- 
Similar solutions also exist for the lower two-component block Hamiltonian, eo(k) + 772 x 2 - straightforward to 
show that the solutions to the eigenvalue problem for the lower block are the same as for the upper one, after one 
makes the replacement kx —>■ —kx- Graphically, these solutions are mirror images of the arcs in Fig. [2] 

Before concluding this subsection, let us also note that the description of the Fermi arc states on the top surface is 
similar to the bottom ones. By assuming that Weyl semimetal is at y < 0 and the vacuum is at y > 0, the appropriate 
boundary conditions are implemented by using the y-dependent parameter m{y) = m9{—y) — fh9{y) and taking the 
limit m —> 00 at the end. Up to a reflection kx —^ — fcx, the corresponding final results for the Fermi arcs on the top 
surface look similar to those on the bottom surface, shown in Fig. [2] 











C. Effective Hamiltonian for surface Fermi arc states 


Following the usual approach in the studies of topological insulator [s^, it may be natural to derive an effective 
Hamiltonian for the surface Fermi arc states. The block Hamiltonians in the simplified model at hand can be naturally 
separated into two parts, i.e., = Hq+H^ , where the zeroth order part Hq corresponds to the original Hamiltonian 

at kx = kz = 0 , i.e., 


TT _ f C'o - C2d^ - 7 m + iVhdl + vdy \ 

° ~ V “" 1 ™ + M2dy - vdy Co - C2dy J ■ 

while Hi contains all the terms with nontrivial dependence on k^ and kz, i.e., 

TT± _ ( Cikl + C2k‘l =F vkx 'jkl - M2k‘l \ 

^ V ~ ^2kx Cikl + C2kl ± vkx ) ' 


(26) 


(27) 


As in the previous analysis, we used ky = —idy. To start with, we have to solve the eigenvalue problem with the 
zeroth order Hamiltonian, i/o'f'o = A^'q. By following the same approach as in Appendix I A H but with kx = kz = 0, 
we find straightforwardly the explicit solutions for the surface Fermi arcs 'Hq. The corresponding energy parameter is 
found to be A = —0.13425 eV. Then, the effective Hamiltonian for the surface states is obtained by integrating over 
the perpendicular direction y, i.e., 


H. 


surf 


foo 1 _ q2 q 

A + y dy't'Ji/i'I'o = A + Cifcj + C2k^ =F vkx q2 + ~ ^2k^) ^ ^ q2 

At “^sur fkx T 7surf^2- 


(28) 


where Q ~ 0.5115, Wsurf ~ 1.440 eV A and 7 surf ~ 17.38 eV A^. Note that the quadratic term in kx vanishes after the 
model parameters are used. 

As is easy to check, the effective Hamiltonian in Eq. (1251) reproduces almost perfectly the shape of the Fermi arcs 
in the kxkz plane. However, it does not contain the information about the finite length of the arcs. We could explain 
this fact in part by pointing out that the corresponding information is encoded in the terms quadratic in momenta 
kx and kz- When such terms are omitted from the zeroth order Hamiltonian Hq, the existence of the surface states 
formally appears to be unconstrained. Therefore, the effective Hamiltonian in Eq. (I28|) will be truly useful only when 
supplemented by its range of validity in the kxkz plane. This, however, seems to diminish its practical value because 
the corresponding range depends on the energy. 


IV. FERMI ARCS IN REALISTIC MODEL 

In this section we will consider the complete low-energy theory described by Hamiltonian m with a 0 . By 
performing a unitary transformation in Eq. o, defined by Uy = ■^^2 ® {I 2 +io'y), we arrive at the following 
equivalent form of the Hamiltonian: 

H = [C2ikl-d^) + Cikl + Co]l2(S)l2-M2{kl-dl)l2®(Jx 

—vkx vdy + ^{k1 — m) —akz{kx — dyY 0 

—vdy T 7 (fcz — rn) vkx 0 akz{kx — dy)^ 

-akz{kx + dyY 0 vkx vdy + ^{kl - m) 

0 akz{kx + dyY —vdy-\-^{k‘l — m) —vkx 

By introducing the spinor wave function T = (^ 1 ,'02, V' 3 j reduce the eigenvalue problem {H — 75)41 = 0 in 

the semimetal (y > 0 ) to the following system of equations: 

[C' 2 (fc^ - dl) - vkx + Cikl T Co - E] 01 T [-M2{kl - dy) + vdy + - 7 m] 02 - akzikx - dy^tpa = 0, (30) 

[-M2ikl - dy) - vdy + -fkl - 7 m] 01 T [C 2 (fc^ - 9^) T vkx + Cikl + Cq - E] 02 T akz{kx - 9 y )^04 = 0, (31) 

-akzikx + dy)'^ipi + [C 2 (fc^ - dy) + vkx + Cikl T Cq - C] 03 T [-M 2 (fc^ - d^) T vdy T ^kl - 7 m] 04 = 0, (32) 

akzikx + 9 y )^02 T [-M 2 (fc^ - d^) - vdy + 'ykl - jm] 03 T [C 2 (fc^ - d^) - vkx + Cikl T Cq - C] 04 = 0. (33) 

On the vacuum side (y < 0), the corresponding set of equations has the same form, but with m replaced by —m. The 
corresponding full set of equations should be also supplemented by the conditions of continuity of the wave functions 
and their derivatives across the vacuum-seinimetal interface at y = 0, see Eqs. (IA16I) through (jA20l) in Appendix I A 21 
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As shown in Appendix lA 2[ the spinor structure of the solution on the semimetal side takes the form: 


2 

«'y>o(2/) 



o-PiV 


(34) 


where the explicit expressions for the exponents are given in Eq. (IA23I) . In the case of surface Fermi arc solutions, 
the exponents take real values. A nontrivial solution exists when the following condition is satisfied: 

{Qt - Qt) {QI - Q 2 ) - {Tt - T}) (Tf - T-) = 0, (35) 

where, by definition, Qf = Q{pi^±kx) and = T{pi,±kx), and the functions Q{p,kx) and T{p,kx) are defined in 
Eqs. (IA12|) and (IA26I) . respectively. 

By taking into account that Tij), kx) vanishes at a = 0, one finds that the above condition reduces to its analog 
in Eq. (1251) in the two-component model. Indeed, a nontrivial solution exists in the model with the two-component 
upper (lower) block Hamiltonian when = Q 2 ) is satisfied. We would like to emphasize that the 

classification of the arc states remains essentially the same also in a general case with a ^ 0. However, because of the 
mixing between the upper and lower block Hamiltonians, the arcs are labeled by the eigenvalues of the operator, 
see Appendix IbI The eigenstates with y = -|-I(y = —I) are the generalizations of the arcs from the upper (lower) 
block Hamiltonian. 

The numerical results for the surface Fermi arc states are shown in Fig. [3] for a = 1 eV A (left panel) and 
a = 50 eV A ( right panel). At fixed energy, there are two surface Fermi arcs related to two different sectors of 
the AaBi (A=Na, K, Rb) compounds with definite eigenvalue of U^. One can check that the wave functions that 
describe these surface Fermi arcs are related to each other by means of the f/H^^ transformation, see Appendix |B] 
By comparing these results with those in the two-component model, see Fig. we find that the quantitative effect 
of a nonzero a on the Fermi arcs is small even when a is moderately large. The only qualitative effect due to a is a 
reconnection of the pair of arcs (from predominantly up and predominantly down sectors) at negative values of the 
Fermi energy. The underlying physics of such an effect is likely to be connected with the loss of the chirality as a 
good quantum number for quasiparticles away from the Dirac/Weyl nodes. Because of the discrete ud-parity, which 
is preserved even at large values of a, there are still two sectors of the theory and there are still small nontrivial arcs 
present, as we see from the right panel of Fig.|31 It will be interesting to explore whether the reconnection of the pairs 
of arcs would also appear in the microscopic theory. It may well be an artifact of the low-energy theory used here. 



FIG. 3: (Color online) The Fermi arcs solutions in the plane of transverse momenta for a = 1 eV A^ (left panel) and a = 50 eV A^ 
(right panel). 
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V. FERMI ARCS AND WEAK BREAKING OF TIME-REVERSAL SYMMETRY 


As we discussed in detail in Sec. Ill B[ the low-energy effective Hamiltonian ([T]) is invariant under the time-reversal 
and inversion symmetries. Moreover, these symmetries play an important role in defining the physical properties of 
A 3 Bi semimetals. Thus, it is natural to ask about possible effects on the structure (and perhaps even the existence) of 
surface Fermi arcs due to breaking of these symmetries. From the physics viewpoint, for example, the corresponding 
discrete symmetries could be broken explicitly by magnetic doping or an external magnetic field. 

In order to study the symmetry breaking effects, we will add to the low-energy Hamiltonian o two additional 
terms controlled by parameters mi and /li: 

= H(k} - ( “ ) ■ ( 36 ) 

' ^ V 0 -fj-ih - (Tzjrni J ' 


By analyzing the Schwinger-Dyson equation for the quasiparticle propagator in AaBi semimetals in a magnetic field, 
we found that these terms are indeed perturbatively generated. Alternatively, these terms can be induced by magnetic 
doping. The value of jii could be interpreted as a mismatch between the chemical potentials of quasiparticle states 
in the Weyl sectors of the theory. The value of mi is a mismatch of the parameter m that determines the chiral shift 
in the two sectors. This means that whenever these symmetry breaking parameters appear, the Z 2 Weyl semimetal 
will get automatically transformed into a true Weyl semimetal with four non-degenerate Weyl nodes. 

By performing a unitary transformation in Eq. (I36L defined by matrix Uy = ® {I 2 + we arrive at the 

following equivalent Hamiltonian: 


Hsh = [C2{kl - dy) + Cikl + Co] I2®l2- M2{kl - dy)l2 ® 

/ -vkx - fli vdy +'^ (kl - m - mi) -akz{kx - dy)^ 

^ —vdy+^[kl — m — mi) vk^ — jj-i 0 

-akz{kx + dyY 0 vkx + Mi 

\ 0 akz(kx + dyY -vdy + ^ [kl - m-\-mi) 


0 

akz{kx - dyY 

vdy + 7 — m -I- mi) 

-vkx + ill 


\ 


(37) 


It is straightforward, although tedious to repeat the same analysis as in Sec. |TVl 

The general surface state solution is of the same type, i.e., ^y^o(y) = , where ff'o = (a, b, c, d)'^ is a constant 

spinor. However, the characteristic equation is considerably more complicated. 


{ [-C2(P^ - kl) + Cikl + Co - Ml - C]' - 
X { [-C 2 (p 2 - kl) + Cikl + Co + Ml - C]" 

+Aa^kl(p^ - kl)^ (m? - y^m?) = 0 . 


[M 2 (p^ - kl) + -/{kl - m - mi)] ^ -h v^{p'^ - kl) - o^kKp^ - klf^ 

- [M2{p^ - kl) + -/{kl -m + mi)Y + v^{p'^ - kl) - a^kl{p^ - A:^)^| 

(38) 


The important effect of the symmetry breaking terms with nonzero mi and mi is that the new characteristic equation 
has four (instead of two degenerate) pairs of distinct solutions: p = ±pi, with i = 1, 2,3,4. The general spinor solution 
in the semimetal takes the following form: 


4 

^!/>o( 2 /) = X! 

i=l 



e-PiV^ 


(39) 


By making use of the equation 






^ _ -C2{pI - kl) + Cikl + Cq - fii - E - vkx ^ _ akz{Pi + kx)^ _ 

—M2{pI — kl) —-/{kl — m — mi) + vpi * — M2{pI — kl) —-/{kl — m — mi) + vp^ 

^ _ _ akz{pi -kx)"^ _^ -C2{pI - kl) + Cikl + Cp + fti - E - vt 

-M2{pI - kl) --/{kl - m + mi) + vpi * -M2{pI - kl) --/{kl - m + mi) + vpi 


(40) 

(41) 


In order to avoid a possible confusion, let us emphasize that the remaining two components Oi and Ci are not 
independent, but fixed unambiguously for each The final solutions for the Fermi arcs are determined after all 
four independent parameters (e.g., with i = 1,2, 3,4) are fixed by satisfying the continuity conditions for the wave 
function at the surface of the semimetal. The corresponding solutions can be obtained by numerical methods. 
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To slightly simplify the analysis, let us consider a special case of vanishing a in more detail. In this case, the states 
from the two-component upper and lower block Hamiltonians decouple. Also, the characteristic equation factorizes, 
effectively giving two separate equations, i.e., 

[-C 2 ip^ - kl) + C,kl + Co - Ml - i?]' - [M 2 (/ - kl) + 7(fc2 - m - mi)] ' + - kl) = Q (up), (42) 

[-C 2 {p^-kl) + Cikl+Co + pi-Ef -[M 2 {p^ -kl)+-i{kl-m + mi)\' + -kl) = 0 (down), (43) 

cf. Eq. (IA8I) . Then, the analysis of the surface Fermi arcs follows very closely the analysis in Sec. IIII 51 






akx/;r 


akx/;r 


akx/;r 


akx/;r 


(a) mi = 10 Ml = 0.05 (b) mi = 10 jl\ = —0.05 (c) mi = 0.005, mi = 10 (d) mi = —0.005, mi = 10 



FIG. 4: The Fermi arcs solutions (thick black lines) in the model with the symmetry breaking parameters mi and fli a.t E — 0. 
The shaded regions represent the projections of the bulk Fermi surfaces onto the plane. The values of mi and jh are 
given in units of A and eV, respectively. 

A number of representative numerical solutions for the Fermi surface arcs in the model with the symmetry breaking 
parameters mi and mi are shown in Fig. |4l The results are obtained for the Fermi energy E = 0. In order to shed 
light on the origin of the individual arcs, in the same figure we also show the projections (shaded regions) of the bulk 
Fermi surfaces onto the k^kz plane. Such a representation reveals that some of the Fermi arcs link disconnected sheets 
of the bulk Fermi surface [l^, while others link different points of the same bulk Fermi surface sheet. 

As suggested by the physical meaning of the symmetry breaking parameters, mi and mi, the Fermi surface arcs for 
the up and down Weyl sectors of the theory are not transformed into each other by a mirror symmetry. In addition 
to the expected effects of (i) changing the length of the arcs (primarily due to nonzero mi) and (ii) shifting the arcs’ 
position in the kx direction (primarily due to nonzero mi), we also see some qualitative changes in the shape and 
branching of the arcs. By comparing Eqs. (HU) and (H51) for the two sectors of the theory, we find that the whole 
asymmetric sets of the Fermi arcs turn into their mirror reflections when both parameters mi and mi change their 
signs. Examples of two pairs of such mirror configurations are shown in panels (e)-(f) and (g)-(h) in Fig. 21 [Strictly 
speaking, the other two pairs of configurations, see (a)-(b) and (c)-(d), are not exact mirror reflections of each other 
because one of the symmetry breaking parameters does not change the sign. Because of a smallness of the parameter, 
there is an appearance of approximate mirror configurations.] 

It is interesting to point out that different topologies of the global (bulk-plus-arcs) Fermi hypersurfaces, including 
the bulk sheets and the surface Fermi arcs, are possible. For example, for a range of symmetry breaking parameters, 
represented by panels (c), (d), (e) and (f) in Fig. 21 we find that the global Fermi hypersurfaces consist of pairs of clearly 
disconnected parts. This is in contrast to the configurations in panels (a) and (b), where different parts touch at four 
points, and in contrast to the configurations in panels (g) and (h), where all parts of the global Fermi hypersurfaces 
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are linked by the Fermi arcs. If samples with completely disconnected parts of the global Fermi hypersurfaces are 
indeed possible, they will be very interesting to study in experiments. 

As we see from panels (g) and (h) in Fig. 01 there are also qualitatively new types of the Fermi arcs possible for a 
range of symmetry breaking parameters. In particular, we find a pair of “short” branches of the Fermi arcs that split 
off from the usual “long” arcs. To the best of our knowledge, the corresponding short arcs have not been predicted 
before. So far, we could not establish a general criterion for the existence of the short arcs. In the configurations in 
panels (g) and (h), they play a profound role by linking two disconnected sheets of the bulk Fermi surface. 


VI. CONCLUSION 

In this paper, we studied the surface Fermi arc states by employing a continuum low-energy effective model. The 
use of analytical methods and a realistic low-energy model provide a deeper insight into the physical properties and 
characteristics of the surface Fermi arcs. In particular, we were able to classify the Fermi arcs with respect to the 
ud-parity and reconfirm the Weyl structure of AaBi semimetals [s^. In this context, it should be noted that the 
experimental observation of the corresponding Fermi arc states have been recently reported for NaaBi [s^. While in 
agreement with the claimed topological semimetal structure, such an observation does not confirm it unambiguously. 
That is because the Fermi arc states are also possible in Dirac materials where the Z 2 Weyl structure is absent [a, [31| . 
The unambiguous confirmation of the I 2 Weyl structure could, however, be established via the quantum oscillations, 
whose period should dependent on the thickness of the semimetal in the same way as in true Weyl semimetals [sil . . 

By introducing the effects of several possible symmetry breaking terms, we show that the Z 2 Weyl structure of AaBi 
is destroyed in a very special way: the compounds become true Weyl semimetals. We suggest that this finding can be 
tested in experiment. For example, by taking into account that the mirror-symmetric pairs of surface Fermi arcs in 
clean AaBi get distorted upon the introduction of explicit symmetry breaking (e.g., by magnetic doping), a number 
of specific features (size, shape and number of branches) should be seen in the surface Fermi arcs. The corresponding 
properties could be studied, for example, by analyzing the quantum oscillations sensitive to the surface states of this 
type M- In the absence of symmetry breaking, there will be a unique period of oscillations dependent in a specific 
way on the thickness of the semimetal slab |32j| . On the other hand, the breaking of symmetry will produce pairs 
of inequivalent arcs of different lengths and the observation of two incommensurate periods of oscillations will be 
expected. In principle, by making use of the analytical results in this study, the details of the oscillations could be 
used to estimate the magnitude of the symmetry breaking terms. 
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Appendix A: Derivations of surface Fermi arcs solutions 

In this Appendix we present the key technical details of deriving the surface Fermi arcs solutions in the 2x2 model, 
introduced in Sec. IIIIBi and in the 4x4 model, introduced in Sec. IIVI 


1. Surface Fermi arcs in 2 x 2 model 

Let us start with the analysis of the surface Fermi arc states in the 2x2 model, introduced in Sec. IIIIBI The 
problem reduces to solving the eigenvalues problem given by Eqs. (|2^ and (l23l) at y > 0 (semimetal), as well as a 
similar set of equations at ?/ < 0 (vacuum), but m replaced by —rh. The corresponding set of equations should be 


also supplemented by the boundary conditions at the vacuum-semimetal interface, i.e., 

^i(-O) = ■*/'i(+0), (Al) 

^2(-0) = M+0), (A2) 

-C2dy'4!i{-0) + M2dy'ii)2{-0) = -C2dyipi{+0) + M2dyilj2{+0), (A3) 

-M29yVji(-0)-I-C'29yV'2(-0) = -M29yV'i(+0) + C'29yV'2(+0), (A4) 
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where ' 01 , 2 ( 2 /) correspond to the vacuum region at y < 0. 

Inside the semimetal (y > 0), the surface state solutions should have the following form: 


^y>o(2/) = ( 6 ) 


(A5) 


By substituting this ansatz in Eqs. 
components a and b: 


and (I23L we arrive at the following set of linear equations for the spinor 


[C 2 {kl - p^) - vkx + Cikl + Cq - E] a + [-M 2 {kl - p^) -vp + - ym] 6 = 0, 

[—M2{k1—p^) + vp + ^k1 — ^m\a+[C2{k1—p^)+vkx+Cik1 + Co — E\b = 0 . 

A nontrivial solution exists when the following characteristic equation is satished: 

[-C'2(/ - kl) + Cikl + Cq-E]^ - [M2{p^ - kl) + 'jkl - + v'^{p‘^ - kl) = 0. 

The solutions to this equation are p = ±pi and p = ±P 2 , where 


Pi = \ kl- 


X + VA2 + Y 


2 (M| - C|) ’ 

Here we introduced the following shorthand notations: 


P2 = \ kl- 


X - VA2 + Y 


2(M| - Cl) 


X = 2C2iCikl + Co-E) + 2jM2{kl-m)-v^, 

Y = A{Ml - Cl) [{C^kl +Co- A)2 - _ 

The spinor components a and 6 of the corresponding nontrivial solution satisfy the constraint 

b -C2{p^ - kl) + Cikl + Co - E - vk, 

a -M 2 {p‘^- kl)--f{kl-m) + vp ' 


(A6) 

(A7) 


(A8) 


(A9) 


(AlO) 

(All) 


(A12) 


Inside the semimetal (y > 0), the wave function should fall off with increasing y. Thus, we use only the negative 
exponents in the general solution, i.e., 




(A13) 


where Qi = Q{pi, kx) with i = 1,2. 

In order to find the vacuum solution (y < 0), we replace m —>■ —m and take the limit m —>■ oo. This leads to the 
following general solution on the vacuum side: 


^y<o(2/) - 





(A14) 


where, for convenience, we took the overall constants to be inversely proportional to The exponents in the 

vacuum solution are determined by 


I 7TO . I 7TO 

It is interesting to note that the conditions of the wave function continuity, see Eqs. (EB) and (IA2I) . are the only 
important conditions to be satisfied. Indeed, the nontrivial solution of Eq. dSB in the limit to —>• oo implies that 
oi = —02 0. This is consistent with Eq. (IA2I1 only when Qi = Q 2 . Concerning the remaining boundary conditions 

in Eqs. (IA3I1 and (IA4I1 . enforcing the continuity of the wave function derivative, they do not add any additional 
constraints. In fact, they are needed only for determining the vacuum spinor components oi and 02 in terms of the 
nontrivial components oi and 02 in the semimetal. Such (finite) solutions always exist. However, as is clear from 
Eq. (IA14I) . the vacuum solution have no much physical content because it vanishes in the limit to —> 00 . 

In conclusion, the boundary conditions at y = 0 are satisfied and, therefore, a nontrivial solution exists when 
Qi = Q 2 - The explicit form of the corresponding condition is given in Eq. (I25D in the main text. 
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2. Surface Fermi arcs in 4 x 4 model 


The analysis of the realistic 4x4 model introduced in Sec. llVl is slightly more involved, but qualitatively similar. 
The eigenvalues problem in this case is given by Eqs. (1501) through (1551) at y > 0 (semimetal), as well as a similar set 
of equations at ?/ < 0 (vacuum), but with m replaced by —fh. The conditions of continuity of the wave functions and 


their derivatives across the vacuum-semimetal surface at y 

^i(-O) = 

-C'29yV'i(-0) -I- M2dytp2{-0) - akzdy'tjjsi-O) = 

M2dy'tjji{-0) - C2dy'tjj2{-0) + akzdy’ijjii-O) = 
-akzdy'^i{-0) - C'29yV'3(-0) -I- A'hdy'ipii-O) = 
ak^dy'4!2{-0) + JVhdy’ijjsi-O) - C2dy'iiji{-0) = 


= 0 are given by 

V'*(+0), for i = l,2,3,4, (A16) 

-C2dyij;i{+0) + M29yV'2(+0) - ak^dyips^+O), (A17) 
M2dyipi{+0) - C2dyilj2i+0) + akzdyijj4{+0), (A18) 
-akzdytpi{+0) - C2dytp3{+0) + M29yV'4(+0), (A19) 

Q;fc^5y'02(+O) -I- M2dy'ijj3{+0) - C2dy'ijj4{+0). (A20) 


In the semimetal (j/ > 0), we look for a general surface state solution in the form 


'^v>o{y) = 


(a 

b 

Vd, 


(A21) 


A nontrivial solution of this type exists when p is a solution to the following characteristic equation; 


[-C'2(p^ - kl) + Cikl + Co-E]^ - [M2(p^ - kl) + -ikl - ym]^ + v'^{p^ - kl) - a^kl{p^ - klf = 0. 


(A22) 


(Strictly speaking, the characteristic equation has the square on the left hand side, implying that the degeneracy of 
its solutions should be doubled.) This equation has two pairs of distinct solutions: p = ±pi and p = ±P 2 , where 


Pi,2 = \ kl- 


x±^/WTY 

2 (M|-C| + a 2 fc 2 )’ 


(A23) 


cf. Eq. (IA9I) . Here, the expression for X is the same as in the 2x2 model (a = 0) in Eq. (lAlOl) . but the expression 
for Y is slightly different, i.e., 


y ^ 4(M2 - + a^kl) [{C4kl + Co- Ef - ^\kl - mf] . 


(A24) 


Therefore, in the half-space occupied by the semimetal (y > 0), the wave function should have the following general 
form: 


^y>o( 2 /) - 


/ ai \ 


/ 02 \ 

Q+ai - T+ci 

g-piy ^ 

Qta2 - T+C2 

Cl 


C2 

\ -yf ai EQiCi y 


\ -T 2 02 + Q^C2 / 


(A25) 


where we also introduced the shorthand notation: = Q{pi, Ek^) and = T(pi, Ek^). Here the function Q(j), kx) 

is the same as in Eq. (IA12I) and 


T{p, kx) = 


ak;,{p + kx)'^ 


—M 2 {p^ — fc2) — 7 (/c 2 — m) + vp 


(A26) 


In the vacuum solution (y < 0), we replace m —>■ —fh and take the limit m —^ oo. In this case, a simple analysis 
leads to the following solution: 


^y<o( 2 /) - 


/ 7 m 


/ ^ “1 \ 
QlC'2+Qfc^Ci 

Cl 

akzQi-\-C2Ci 

V yJCl-cflkl } 


pPiy 


-k 




( a2 \ 

_ Q 2 C* 2 +Q!fczC 2 
^Cl-a'^kl 
C2 

Q:fc^Q.2+C*2 C 2 

V y/ci-a'^kl } 


^P2y 


(A27) 
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where, for convenience, we introduced an overall constant inversely proportional to \Jym. The explicit form of the 
exponents in this solution is determined by 


Pi ^ 


'ym 


-M 2 + v'C'l - a2fe2 ^ 


P2 ^ 


7 m 


-M 2 - 


(A28) 


Note that the signs in the exponents of the vacuum solution (IA27I1 are chosen so that the wave function vanishes at 
y —>■ — 00 . 

The conditions of the continuity of the wave function in Eq. (IA16I) lead to the following constraints: 

0 = Qfai - T^ci + Qta2 - T+C 2 , (A29) 

0 = Tj“fli — Cl + T 2 a2 — Q 2 C 2 , (A30) 

together with 02 = —oi and C 2 = —ci. Here, we took into account that the left hand side of Eq. (jAl 6 () vanishes in 

the limit rh ^ 00. 

As in the case of a two-component model, discussed in Appendix [xn there is no need to satisfy the continuity 
conditions for the wave function derivatives, given by Eqs. (IA17I) through (IA20I) . The reason is that these conditions 
add no additional constraints on the spinor solutions in the semimetal. They are needed only for determining the 
components of the vacuum solution at y < 0. Since the latter has no physical content in the limit to —> 00 , we can 
safely ignore the conditions in Eqs. (IA17I) through (IA20I) . 

In order to have a nontrivial solution to Eqs. (IA29I) and (IA30I) . the condition in Eq. (1351) in the main text should 
be satisfied. 


Appendix B: Symmetries and surface Fermi arcs bispinors 


In this Appendix, we discuss the properties of the Fermi arc surface states with respect to the discrete symmetries 
= C/Hfc^ and , introduced in Sec. Ill B1 To start with, let us note that the general Fermi arc spinor in Eq. (IA25I) 
contains all possible solutions. It is possible to classify these solutions with respect to the discrete symmetry by 
choosing them as eigenstates of the operator U^. 

In order to construct the first group of solutions, we use the relations in Eqs. (IA29() and (jA30() and rewrite the 
Fermi arc spinor in Eq. (IA25I) in the following form: 


T+ = ai 



( ^ ^ 


( 1 \ 


P+ 7^+ "^1 ^2 

- ai 

rp-\- 

ai 

T -TP 

1 

1 


Q\ “Q2 


Qi Q 2 


T~ + 0~ 

V Qr-Q^"/ 


T~ + 0~ 

V 2+b/2Q^_Q-y 


(Bl) 


By noting that Pi’s (with i = 1,2), defined in Eq. (IA23I) . contain only quadratic terms in momenta and we 
conclude that both of them are invariant under the H^^ and Hfc^ transformations. The other quantities, used in 
Eq. (IBII) . transform as follows: 


Bfc, Qf = Qt , Hfe^ = -7;±, 


(B2) 

(B3) 


It is straightforward to check that the spinor in Eq. (IBID is an eigenstate of the operator with the eigenvalue 
X = +t- Indeed, by making use of the definition of the matrix U, we find that 


/ 


U'i’+ = ai 


V 


Qt - T{ 


-t: 


1 


Q 


TP -TP 


^ Qi Q 2 


\TP-Tp^ 


e-P^y - ai 


\TP-Tpl 

Qi —Q 2 . 




Qi —Q 2 . 




TP-TP 


-t: 


Q 


- TP -TP 


^ Qi Q 2 . 


o-P^y 


(B4) 


Then, by taking into account that we see that ~ claimed. 
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By using the relations in Eqs. (IA29|) and (jA30l) . the Fermi arc spinor in Eq. (IA25I) can be also rewritten in the 
following alternative form: 


/ 


= Cl 


77-77 

Qt-Qi 

'r- 

1 


\ 


+ 


/ 


Cl 


V "^1 y 


T,+ -Td 


\ 


Qt-Qt 

/0+ T^ + 

1 

V ■^2' + Q2 ) 


(B5) 


In this case, as is easy to check, the spinor is an eigenstate of the operator with eigenvalue y = — 1. Indeed, by 
making use of the definition of the matrix U, we find that 


/ 


[/«'- = Cl 


rp-f- rp-\- 
-^1 -^2 

Qt-Qt 


Ql 


+ A 


Qt-Qt 


-T, 


-1 

1 Qt-Qt 


Ql 


\ 


( 


o-PlV . 


Cl 


Ty-Ty 


Qt-Qt 

fi+'^t-'^'t 7 ^+ 

^2 


-T. 
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which implies that In other words, the eigenstate corresponds to x = —1, as claimed. 

Now, let us explore the implications of the UWk^ symmetry in the model at hand. By applying the corresponding 
operator to the eigenstates 'I'i, we arrive at the following results: 
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These results show that 4'± are not eigenstates of the operator UHk^. However, by taking into account the constraint 
in Eq. (1551) . one can check that the operator UHk^ interchanges the two types of the states, i.e., 4'+ 

In conclusion, the results of this Appendix confirm the claim in the main text of the paper concerning the symmetry 
properties of the low-energy theory for AsBi (A=Na, K, Rb) semimetals, as well as the classification of their Fermi arc 
states. These are in complete agreement with the claim that the corresponding compounds are Z 2 Weyl semimetals. 
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